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Abstract 

We investigate the multi-photon process of above-threshold ionization for the light elements hy¬ 
drogen, carbon, nitrogen and oxygen in the hard x-ray regime. Numerical challenges are discussed 
and by comparing Hartree-Fock-Slater calculations to conhguration-interaction-singles results we 
justify the mean-held potential approach in this regime. We present a theoretical prediction of 
two-photon above-threshold-ionization cross sections for the mentioned elements. Moreover, we 
study how the importance of above-threshold ionization varies with intensity. We hnd that for car¬ 
bon, at x-ray intensities around lO^^Wcm”^, two-photon above-threshold ionization of the K-shell 
electrons is as probable as one-photon ionization of the L-shell electrons. 
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I. INTRODUCTION 


The development of x-ray free-electron lasers (XFELs) in recent years has allowed the 
production of ultrashort x-ray pulses at ever increasing intensities. Highly intense hard x 
rays are of particular interest, e.g., for the purposes of molecular imaging at atomic resolution 
[1]. Anticipating further developments in the direction of ultrashort x-ray pulses down to a 
few hundred attoseconds time scale (1 as= 10“^® s), XFELs will represent the ideal tool for 
single-molecule imaging via coherent x-ray scattering Further interesting applications 

of highly intense and ultrashort pulses include the investigation of electronic dynamics in 
atoms and (bio-) molecules which typically take place on a time scale between attoseconds 
and tens of femtoseconds s) [11[5]. 

The Linac Coherent Light Source (LCLS) XFEL has been running since 2009 at the SLAG 
National Accelerator Laboratory in the US and was the first XFEL capable of producing hard 
x-rays [6]. Recently, the SACLA XFEL at the SPring-8 facility in Japan reached intensities 
of 10^° Wcm“^ at photon energy of 9.9 keV [7]. The European XFEL in Hamburg, Germany, 
is intended to produce x-rays with photon energies up to 12 — 100 keV [S]. By focusing the 
pulses down to a few nanometers, XFELs can reach intensities that are orders of magnitudes 
greater than previously achieved. Even though the interaction probability of x-rays with 
matter is low |9], in this high-intensity regime it is necessary to consider the importance of 
nonlinear processes affecting electronic dynamics of atomic, molecular or solid-state target 
systems. 

The present study focuses on the nonlinear effect of above-threshold ionization (ATI). 
First observed in 1997 by P. Agostini et ai, ATI is a process whereby an electron absorbs 
more photons than are necessary for ionization [10]. This process has been studied exten¬ 
sively by now, especially in the range of infrared to visible and XUV light [nHl3]; however, 
less work exists for x-ray ATI in the high-intensity regime [T3|, and mostly on hydrogen or 
hydrogen-like ions [T5IIT8] . 

The purpose of this work is to examine the role and the magnitude of ATI in the x-ray 
regime under the high-intensity conditions that will become available soon at XFELs. To this 
end, photoelectron spectra are calculated in order to quantify the effect of ATI. Our method 
for the calculation of photoelectron spectra is based on the time-dependent configuration 
interaction singles (TDCIS) method [T^IET] . The A^-electron wave function is obtained by 
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solving the Schrodinger equation numerically on a grid. In order to prevent reflections from 
the end of the grid and to analyze the outgoing wave packet, the so-called wave function 
splitting method is employed. The first-principle calculations of the TDCIS method are 
compared to calculations treating the atomic potential on the Hartree-Fock-Slater (HFS) 
level. For the purposes of molecular imaging the most important elements to consider are 
those commonly found in organic molecules. For this reason carbon, nitrogen, oxygen and 
hydrogen were chosen as the focus of this study. The two-photon ATI cross sections for 
these elements were calculated at the representative and commonly used hard x-ray photon 
energies of 8, 10 and 12 keV. An investigation was also carried out to see at which intensities 
ATI makes an important contribution to the overall ionization probability. 

In Sec. [IT] a discussion of the theoretical basis of the method is presented, followed by 
the numerical and computational challenges faced during the investigation of photon-atom 


interactions in the x-ray regime. In Sec. Ill both the two-photon ATI cross sections for 
light elements and the results of an intensity study for ATI in carbon are shown. A short 


summary in Sec. IV concludes this work. Atomic units are used throughout unless otherwise 
stated. 


II. THEORY AND METHOD 

A. Theoretical Overview 

The time-dependent Schrodinger equation for a A-electron system is given by 

!k|»(«)) = F|*(«)>. (1) 

Here |T(f)) is the full Y-electron wave function. The Hamiltonian H{t) is of the form 

H{t) = Ho + Hi + p-A{t), ( 2 ) 

where Hq = T + V^uc + Vmf — Ehf, containing terms for the kinetic energy T, the nuclear 
potential Kiuc, the mean-held potential IVif, and the Hartree-Fock energy F^hf- The part 
of the exact electron-electron Coulomb interaction V not included within is provided 
by iFi = 14 ~ fmf- Under the dipole approximation p ■ A{t) describes the light-matter 
interaction in the velocity form (the light is assumed to be linearly polarized). 
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The solution to Eq. Q is written as a superposition of the Hartree-Fock ground state 
|<ho) and all 1-particle-l-hole excitations |<h“) 

|>I>(())=a„M|4„) + 5^<M|4>“). (3) 

i,a 

where the index i specifies an initially occupied orbital and a indicates an unoccupied orbital 
to which the electron can be excited. The 1-particle-l-hole excitations are given by 

l^“) = ^ {ci+Ci+ + cl_Ci_| \%). (4) 

Here, and c represent creation and annihilation operators for the corresponding orbital, 
respectively. By applying the time-dependent Schrodinger equation to Eq. ([^ and projecting 
onto the |<ho) and |<1>“) states we obtain an equation of motion for the expansion coefficients 
ao{t) and a“(t). The TDCIS method as described here assumes closed shell atoms in the 
Hartree-Fock ground state, such that the total spin is zero. A detailed description of TDCIS 
can be found in Refs. [HI [22] ■ 

In order to both eliminate reflections and to store information about the outgoing wave 
packet, the method of wave-function splitting was used, which was first implemented by 
Tong et al. [23]. Briefly, a splitting operator S is applied to the wave function that has the 
shape of a smoothed-out step function 

S = {l-Fexp[-(f-rc)/A]}”\ (5) 


where the parameter A controls the smoothness of the function and determines the center. 
The wave function is then split into |Xin(to)) and |Xout(^o)) where the inner and outer parts 
are treated separately: 

IXin(^o)) = (1 - S)IXi(fo)), 


( 6 ) 


IXout(io)) = S\Xi{to)). 

Here, to is the current time step and |Xi(^o)) is the wave function for a particular ionization 
channel i. By absorbing and independently propagating the outer part analytically to long 
times, reflections from the end of the grid are avoided. A more detailed description of 
the calculation of photoelectron spectra (PES) using the TDCIS approach can be found in 
Ref. [20]. In order to compute cross sections for open shell atoms, we use here the HFS 
atomic potential [23]. The HFS potentials for the various elements were calculated using 
the XATOM code [25] . 
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B. Numerical Challenges 


During the investigation of the interaction of x-rays with atoms a variety of computational 
and numerical challenges arise. The grid on which the wave function is represented requires 
a large radius in order to efficiently apply the splitting method to the high-energy wave 
packets produced. Grid sizes of around 120 Bohr radii were found to be large enough. The 
number of grid points was chosen at approximately 10 points per de-Broglie wavelength in 
order to well represent the high-energy parts of the outgoing wave function. Increasing the 
number of grid points further did not influence the PES. It was found that for the x-ray 
photon-energy regime choosing a maximum angular momentum of higher than 3 did not 
signihcantly affect the PES. 

Because of the high photon energy the propagation time step used to propagate the wave 
function needs to be very small. We found that using Runge Kutta to the 4th order about 
20 time steps per electric held oscillation are required to prevent signihcant artifacts from 
appearing in the PES. As a consequence of the small propagation time step the splitting 
function had to be applied very frequently. Applying the splitting every 3 time steps is 
appropriate to remove artifacts due to rehections. However, this causes some practical 
problems as large amounts of data must be stored. Certain parameters are found to be 
mostly unimportant for convergence. Consistently with the study in Ref. [20], the smoothing 
of the splitting function A has little to no effect as long as the wave function close to the 
nucleus is not disturbed. 

Obviously, using the HFS potential signihcantly reduces the computational time as it 
removes the necessity to calculate the large number of Coulomb matrix elements, i.e., the 
exact Coulomb interaction between the electrons. In order to hnd out if the HFS approach is 
valid in the x-ray photon-energy range, neon was studied with both methods before applying 
it to the elements carbon, nitrogen and oxygen, which are more relevant for chemistry and 
biology than nobel-gas atoms. In Fig. we present the photoelectron spectrum of neon 
for a Gaussian pulse of 10 keV photon energy and 0.48 fs pulse duration at an intensity of 
10^^ Wcm“^ using the HFS potential (green) and the full TDCIS method (red). Except for 
the slight upward energy shift for HFS the shape and height of the peaks are the same. The 
noise levels and artifacts present in the spectrum are also unaffected by the method used. 
The comparison suggests that in the photon energy regime of hard x-rays the HFS approach 
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Figure 1: (Color online) Photoelectron spectrum of neon for a Gaussian pulse with 10 keV photon 
energy, at 10^^ Wcm“^ intensity and 0.48 fs pulse duration, using HFS (green) and TDCIS (red). 
The one-photon ionization peak the first two ATI peaks are shown. Each peak consists of two 
subpeaks, the one at lower energy being associated with ionization from the K shell, the other one 
corresponding to L-shell ionization. 

is well justified for open shell atoms with a similar Z value, namely carbon, nitrogen, and 
oxygen. Our calculations for neon show that electron correlation effects play a minor role 
for one- and two-photon absorption in this regime. 


III. RESULTS AND DISCUSSION 

We calculated two-photon ATI cross sections for hydrogen and compared them with 
previous work |T3]. They are found to be in very close agreement. Using 10^° Wcm“^, a 
relative difference of 2% was found at the photon energy of 10 keV and 7% at 8 keV. With 
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a slightly larger intensity of 3.5 x 10^° Wcm ^ the cross section at 5 keV was found to be 
2.02 X 10 cm'^s, which amounts to a difference of 1% compared to Ref. [T^. Since the 
method employed here relies on the numerically exact solution of the Schrodinger equation 
Q, no difficulties arise with sums over intermediate states, which appear in a perturbative 
treatment na. Reference [I5] demonstrates that the interaction, which is not included 
here, may be neglected in the photon-energy range of current interest. As already indicated 
in Eq. ([^ by the exclusive time dependence of the vector potential all calculations were 
performed under the assumption of the dipole approximation. However, under the conditions 
of short-wavelength x-rays this assumption may no longer be valid. Indeed, Zhou and Chu 
indicate that nondipole effects signihcantly change the photoelectron angular distribution 
and that the nondipole ATI spectra are enhanced in the high photon-energy regime [26] . 
On the other hand, it was shown for the x-ray regime that when including all multipoles 
the total two-photon ATI cross section differs less than an order of magnitude from the 
cross section calculated in dipole approximation for sufficiently small nuclear charge [TH] . 
Therefore, although our values may underestimate the two-photon ATI cross section, we 
expect this underestimation in the integrated spectrum to be much smaller than an order 
of magnitude (a factor of around 2 — 3) [27] . 

The two-photon ATI cross sections were found in the perturbative limit, i.e., by assur¬ 
ing that the ionization probability is low enough to not deplete the ground state. In the 
perturbative limit the ionization probability due to 2-photon absorption is given by 
^(2)^(2) _ ^( 2 ) j j2^^^ where is the two-photon cross section and J is the photon flux 
in cm“^s“^; is the fluence for two-photon absorption. Assuming a Gaussian pulse the 
electric held has the form 

E{t) = Eo cos(a;f)e-2^“2 (7) 


where E^ is the peak electric held, u is the central held frequency, and r is the pulse duration. 
Then, the huence E^‘^1 is given by 


Fr'4 

f(2) = TA 
a;2 


TT 


( 8 ) 


dn 2 V Stt 

The two-photon ATI cross sections for hydrogen, carbon, nitrogen, oxygen, and neon at 
the hard x-ray energies of 8, 10, and 12 keV are presented in Table |T| We hnd the expected 
increase in the two-photon ATI cross section for higher Z values and a drop with larger 
photon energy. As mentioned previously, all cross sections were found under the dipole 
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Table I: Two-photon ATI cross sections for the light elements (E is the photon energy). All were 
calculated from integrating the corresponding photoelectron peaks at 10^® Wcm“^ intensity and 
0.12 fs of pulse duration. 


E (keV) 

Two-photon cross sections (cm‘^s) 

Hydrogen 

Carbon 

Nitrogen 

Oxygen 

Neon 

8 

1.44 X 10“®® 

1.64 X 10“®’^ 

3.21 X 10“®® 

5.62 X 10“®® 

1.65 X 10“®’ 

10 

4.69 X 10“®’' 

4.61 X 10“®® 

9.23 X 10“®® 

1.70 X 10“®® 

3.71 X 10“®® 

12 

1.72 X 10“®’' 

1.79 X 10“®® 

3.82 X 10“®® 

6.94 X 10“®® 

1.95 X 10“®® 


approximation, which is expected to slightly underestimate the cross sections. The two- 
photon ATI cross section of beryllium was also calculated at 10 keV photon energy to be 
1.25 X 10-®3 cm^s. 

Because imaging of organic molecules is of particular interest we perform an intensity 
study on carbon with an incoming photon energy of 10 keV in order to hnd the regime in 
which the ionization due to two-photon ATI is of the same order as that for one-photon 
ionization. As seen in Fig. at an intensity of 10^^ Wcm“^ the depopulation due to K-shell 
ATI in the direction 6 = 0 becomes higher than valence one-photon ionization. The PES 
in the directions vr/b, vr/S, and 7r/2 are also shown (for better visualization only 15 orders 
of magnitude are shown). Note, that a small peak around zero kinetic energy is observed 
[281 |29]. However, the height of the peak might be underestimated, because due to our 
splitting approach the propagation time must be sufficiently large in order to detect all 
electrons of interest in the splitting region. The results shown in Fig. |^were produced in a 
calculation where the propagation time was approximately 20 times shorter than necessary 
to collect all electrons with a kinetic energy on the scale of 1 eV. In order to elucidate 
this slow-electron peak further, a new calculation was performed, now involving a pulse at 
10 keV photon energy, of 12 as duration and 3.5 ■ 10^^ Wcm“^ intensity in order to be able to 
propagate long enough and to observe the slow-electron peak in its full height. The results 
are presented in Fig. The emergence of the peak can be attributed to the bandwidth of 
the pulse which spans the binding energy of the valence electrons: after the absorption of 
one photon by the valence shell of the atom the emission of a photon can occur and, thereby, 
slow electrons are produced. 

In Fig. [^the ionization probability of carbon is shown as a function of intensity for the 
case of one-photon ionization out of the L shell together with two-photon ionization out of 

















Figure 2; (Color online) Photoelectron spectra of carbon for different angles: a) 0, b) 7r/6, c) vr/S, 
d) 7r/2. The spectra are calcnlated for a pulse centered at 10 keV photon energy, at an intensity 
of 10^^ Wcm“^ and 0.12 fs pulse duration. Note that the height of the K-shell ATI peak is greater 
than or comparable to the one-photon L-shell ionization peak. 

the K shell. We see the characteristic quadratic behavior of the ATI peak as a function 
of the intensity and the linear behavior in the one-photon valence ionization probability. 
Saturation effects do not play a role until the intensity range near Wcm“^ is reached. 
In fact, one can see that at an intensity between 10^^ — 10^^ Wcm“^ the fully angle- and 
energy-integrated K-shell ATI peak, i.e., the ionization probability due to ATI out of the K 
shell, is comparable to the probability to ionize with one photon out of the valence shells. 
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Figure 3: (Color online) Depopulation at different intensities due to two-photon ATI against one- 
photon valence electron ionization for carbon, at 10 keV photon energy and 0.48 fs pulse duration. 
The ionization probability was found by integrating the corresponding PES peaks. 

IV. CONCLUSION 

We have presented a prediction for the two-photon ATI cross sections of the light elements 
carbon, nitrogen, and oxygen at hard x-ray energies common for current experiments at 
XFELs. At intensities that one may reach for future hard x-ray experiments, scientists 
should consider how ATI will affect their results. We conclude that ATI remains a negligible 
fraction of ionization for intensities at the most recent XFEL experiments with hard x-rays. 
However, we predict that with photon energies at around 10 keV, when entering the regime 
of 10^^ Wcm“^ and above, the ionization probability of the core electrons by ATI approaches 
the same order of magnitude as valence stripping by one-photon ionization for elements with 
a similar nuclear charge Z as carbon. It is likely that the neglected nondipole effects enhance 
the ATI spectrum by a factor smaller than an order of magnitude. Therefore, we can present 
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Figure 4: (Color online) Photoelectron yield for carbon, shown for a pulse of 12 as duration and 
3.5 • 10^^ Wcni“^ intensity with a photon energy of 10 keV after a sufficiently long propagation 
time. The slow-electron peak (magnified in the inset) as well as the 1-photon peak and the first 
two ATI peaks are shown for 4 different angles with respect to the polarization direction. The 
height of the slow-electron peak is comparable to the first ATI peak in the directions of vr/S and 
7r/2. 

our results as a lower limit on the importance of ATI and suggest that ATI be taken into 
account when entering this high-intensity regime. In particular, we hope that the data 
presented can be a guide for future experiments investigating imaging and nonlinear x-ray 
optics. 
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